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Abstract. We consider a version of a Glauber dynamics for ap-spin Sherrington- 
Kirkpatrick model of a spin glass that can be seen as a time change of simple 
random walk on the iV-dimensional hypercube. We show that, for any p > 3 and 
any inverse temperature /3 > 0, there exist constants 70 > 0, such that for all 
. exponential time scales, exp(7JV), with 7 < 70, the properly rescaled clock process 

(time-change process), converges to an a-stable subordinator where a = j/f3 2 < 1. 
Moreover, the dynamics exhibits aging at these time scales with time-time correla- 
tion function converging to the arcsine law of this a-stable subordinator. In other 



(N 



words, up to rescaling, on these time scales (that are shorter than the equilibration 
time of the system), the dynamics of p-spin models ages in the same way as the 
REM, and by extension Bouchaud's REM-like trap model, confirming the latter 
as a universal aging mechanism for a wide range of systems. The SK model (the 
case p = 2) seems to belong to a different universality class. 

5?" 

^3 ■ I. Introduction and results 

-t— > 

a . 

Aging has become one of the main paradigms to describe the long-time behav- 
ior of complex and/or disordered systems. Systems that have strongly motivated 
this research are spin glasses, where aging was first observed experimentally in the 
■ anomalous relaxation patterns of the magnetization [LSNB83, Cha84j. The theo- 

retical modeling of aging phenomena took a major leap with the introduction of 
so-called trap models by Bouchaud and Dean in the early 1990'ies |Bou92l [BD95J 
(see |BCKM 98j for a review). These models reproduce the characteristic power law 
\Q ■ behavior seen experimentally while being sufficiently simple to allow for detailed 

analytical treatment. While trap models are heuristically motivated to capture the 
behavior of the dynamics of spin glass models, there is no clear theoretical, let alone 
mathematical derivation of these from an underlying spin-glass dynamics. The first 
attempt to establish such a connection was made in [BBG02, BBG03a, BBG03bJ 
where it was shown that starting from a particular Glauber dynamics of the Ran- 
dom Energy Model (REM), at low temperatures and at the time scale slightly shorter 
than the equilibration time of the dynamics, the aging of the time-time correlation 
function of the dynamics converged to that given by Bouchaud's REM-like trap 
model. 



On the other hand, in a series of papers [BC051 IBCM061 lBC07al lBC07bj a sys- 
tematic investigation of a variety of trap models was initiated. In this process, it 
emerged that there appears to be an almost universal aging mechanism based on 
a-stable subordinators that governs aging in most of the trap models. It was also 
shown that the same feature holds for the dynamics of the REM at shorter time 
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scales than those considered in [BBGOSa, IBBG03b] . and that this also hap pens at 



high temperature provided appro priate t ime scales are considered [BC07aj. For a 



general review on trap models see [BC06 



In all models considered so far, however, the random variables describing the 
quenched disorder were considered to be independent, be it in the REM or in the 
trap models. Aging in correlated spin glass models was investigated rigorously only 
in some cases of spherical SK models and at very short time scales |BDG01j . In the 
present paper we show for the first time that the same type of aging mechanism is 
relevant also in correlated spin glasses, at least on time scales that are short compared 
to equilibration time (but exponentially large in the volume of the system). 

Let us first describe the class of models we are considering. Our state spaces will 
be the iV-dimensional hypercube, <Sjv = { — 1, 1}^- Rn '■ <Sn x <Sn —> [— 1) 1] denotes 
as usual the normalized overlap, Rn(o~, t) = N^ 1 YliLi % T i- The Hamiltonian of the 
p-spin SK-model is defined as ^/NH^, where : Sn — > R is the centered normal 
process indexed by with covariance 

E[H N (a)H N (r)] = R N (a,ry, (1.1) 

and p G N, p > 2. We will denote by 7i the a- algebra generated by the random 
variables Hn(o~), a G Sn, N G N. The corresponding Gibbs measure is then given by 

^(ajE^'W, (1.2) 

where Zp^ denotes the normalizing partition function. 

We define the classical trap-model dynamics as a nearest neighbor continuous time 
Markov chain erjv(-) on with transition rates 

, s ( N -l e -PVNH N (a) if dist(d,r) = 1, 

w n (<t,t) = < . (1.3) 

I U, otherwise; 

here dist(-, •) is the graph distance on the hypercube, 

1 N 

dist(cr, r) = ^ \(Tj - Tj\ . (1.4) 

i=X 

A simple way to construct this dynamics is as a time change of a simple random 
walk on S^'- We denote by Y N {k) G <Sjv, k G N, the simple unbiased random walk 
(SRW) on S N started at some fixed point of Sn, say at {1, . . . , 1}. For f3 > we 
define the clock-process by 

fc-i 

S N {k) = J2 e * ex P {^ h n {Y N {i))}, (1.5) 

i=0 

where {ei,i G N} is a sequence of mean-one i.i.d. exponential random variables. We 
denote by y the cx-algebra generated by the SRW random variables Y^(k), k G N, 
N G N. The cr-algebra generated by the random variables e^z G N will be denoted 
by S. Then the process crjv(-) can be written as 

a N (t) = Y N (S^(t)). (1.6) 

Obviously, o~n is reversible with respect to the measure /x^at- We will consider all 
random processes to be defined on an abstract probability space (Q, JF, P). Note that 
the three a-algebras TC, y, and £ are all independent under P. 
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We will systematically use the definition of the dynamics given by (11.31) or (\l.Q\i . 
This is the same as was used in the analysis of the REM and in most work on trap 
models. It differs substantially from more popular dynamics such as the Metropolis 
or the heat-bath algorithm. The main difference is that in these dynamics the tra- 
jectories are not independent of the environment and are biased against going up in 
energy. This may have a substantial effect on the dynamics, and we do not know 
whether our results will apply also (with some modifications) in these cases. The 
fact is that we currently do not have the tools to analyze these dynamics even in the 
case of the REM! 

Let V a (t) be the a-stable subordinator with the Laplace transform given by 

E [ e -AVa(t)] = exp (_ tA «). (1.7) 

The main technical result on the dynamics will be the following theorem that provides 
the asymptotic behavior of the clock process. 

Theorem 1.1. There exists a function ((p) such that for all p > 3 and 7 satisfying 

0<7<min(/3 2 ,C(p)/3), (1.8) 

under the conditional distribution P[-|3^] the law of the stochastic process 

S N (t) = e^ N S N {[tN 1 / 2 e N '' 2 ^ 2 \), t > 0, (1.9) 

defined on the the space of cadlag functions equipped with the Skorokhod M\-topology, 
converges, y-a.s., to the law of '7/ '(3 2 - stable subordinator V^jp(Kt\t > 0, where K 
is a positive constant depending on 7, j3 and p. 

Moreover, the function ((p) is increasing and it satisfies 



C(3)~ 1.0291 and lim ((p) = ^2 log 2. (1.10) 

p— *oo 

We will explain in Section [5] what the Mi-topology is. Roughly, it is a weak 
topology that does not convey much information at the jumps of the limiting process: 
it can be the case that the approximating processes jumps several times at rather 
short distances to produce one bigger jump of the limit process. This will actually 
be the case in our models for p < 00, while it is not the case in the REM. Therefore 
we cannot replace the Mi topology with the stronger Jx-topology in Theorem 11.11 

To control the behavior of spin-spin correlation functions that are commonly used 
to characterize aging, we need to know more on how these jumps occur at finite N. 
What we will show, is that if we the slightly coarse-grain the process Sn over blocks 
of size o(N), the rescaled process does converge in the Ji-topology. What this says, 
is that the jumps of the limiting process are compounded by smaller jumps that 
are made over < o(N) steps of the SRW. In other words, the jumps of the limiting 
process come from waiting times accumulated in one slightly extended trap, and 
during this entire time only a negligible fraction of the spins are flipped. That will 
imply the following aging result. 

Theorem 1.2. Let A £ N (t,s) be the event defined by 

A e N (t,s) = {R N (a N (tei N ),a N ((t + sy N )) >l-e}. (1.11) 

Then, under the hypothesis of Theorem ! 1. 11 for all e G (0, 1), t > and s > 0, 

sin mi 

lim F[A £ N {t,s)} = —— u a -\l-u)- a du. (1.12) 

N^oo IX J 
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Remark. We will in fact prove the stronger statement that aging in the above sense 
occurs along almost every random walk trajectory, that is 

sin nir 

lim F[A E N (t,s)\y] = / u a -\l-u)~ a du, y-a.s. (1.13) 

JV^oo 7T Jq 

Let us discuss the meaning of these results. e lN is the time-scale at which we want 
to observe the process. According to Theorem II. 1[ at this time the random walk 
will make of the order of N 1 ^ 2 e N " / / 2/3 <C e yN steps. Since this number is also much 
smaller than 2 N (as follows from (11.101) ). the random walk will essentially visit that 
number of sites. 

If the random process Hn was i.i.d., then the maximum of if at along the trajectory 
would be {2\n{N 1 / 2 e N ^ 2 l2l32 )) 1/2 ~ iV 1/2 7//3, and the time spent in that site would 
be o f order e yN . Since Theorem 11.11 holds also in the i.i.d. case, that is in the REM 
(see BC07aj ). the time spent in the maximum is comparable to the total time and 
the convergence to the a-stable subordinator implies that the total accumulated time 
is composed of pieces of order e^ N that are collected along the trajectory. In fact, 
each jump of the subordinator corresponds to one visit to a site that has waiting 
times of that order. In a common metaphor, the sites are referred to as traps and 
the mean waiting times as their depths. 

The theorem in the general case states that in the p-spin model, the same is 
essentially true. The difference will be that the traps here will not consist of a 
single site, but consist of a deep valley (along the trajectory) whose bottom that has 
approximately the same energy as in the i.i.d. case and whose shape and width we 
will be able to describe quite precisely. Remarkably, the number of sites contributing 
significantly to the residence time in the valley is essentially finite, and different 
valleys are statistically independent. 

The fact that traps are finite may appear quite surprising to those familiar with 
the statics of p-spin models. From the results there (see [Tal031 IBov06j ). it is known 
that the Gibbs measure concentrates on "lumps" whose diameter is of order Ne p , 
with e p > 0. The mystery is however solved easily: the process H N (a) does indeed 
decreases essentially linearly with speed N~ l l 2 from a local maximum. Thus, the 
residence times in such sites decrease geometrically, so that the contributions of a 
neighborhood of size K of a local maximum amounts to a fraction of (1 — c~ K ) of 
the total time spend in that valley ; for the support of the Gibbs measure, one needs 
however to take into account the entropy, that is that the volumes of the balls of 
radius r increases like N r . For the dynamics, at least at our time-scales, this is, 
however, irrelevant, since the SRW leaves a local minimum essentially ballistically. 

The proof of Theorem 1 1 . 1 1 relies on the combination of detailed information on the 
properties of simple random walk on the hypercube, which is provided in Section H] 
(but see also [Mat89l IBG061 ICG06j ). and comparison of the process on the 
trajectory of the SRW to a simpler Gaussian process using interpolation techniques 
a la Slepian, familiar from extreme value theory of Gaussian processes. 

Let us explain this in more detail. On the time scales we are considering, the SRW 
makes tiV 1/2 exp(iV7 2 /2 / 3 2 ) < tiV 1/2 exp(A^C(p) 2 /2) < 2 N steps. In this regime the 
SRW is extremely "transient", in the sense that (i) starting from a given point x, for 
a times t < v ~ N u , u < 1, the distance from x grows essentially linearly with speed 
one, that is there are no backtrackings with high probability; (ii) the SRW will never 
return to a neighborhood of size v of the starting point x, with high probability. The 
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upshot is that we can think of the trajectory of the SRW essentially as of a straight 
line. 

Next we consider the Gaussian process restricted to the SRW trajectory. We expect 
that the main contributions to the sums Sjv(A;) come from places where Yjv is maximal 
(on the trajectory). We expect that the distribution of these extremes do not feel the 
correlation between points farther than v apart. On the other hand, for points closer 
than v, the correlation function R^{Y]sr{i),Y^{j)) p can be well approximated by a 
linear function 1 — 2p\i — j\/N (using that R^(Y N (i), Y N (j)) ~ 1 — 2\i — j\/N). This is 
convenient since this process has an explicit representation in terms of i.i.d. random 
variables that allow for explicit computations (in fact, this is one of the famous 
Slepian processes for which the extremal distribution can be computed explicitly 
|Sle611 IShe71j ). Thus the idea is to cut the SRW trajectory into blocks of length 
v and to replace the original process H^{Y^{ij) by a new one Ui, where Ui and Uj 
are independent, if i, j are not in the same block, and E[[/j[/j] = 1 — 2p\i — j\/N 
if they are. For the new process, Theorem 11.11 is relatively straightforward. The 
main step is the computation of Laplace transforms in Section [2j Comparing the 
real process with the auxiliary one is the bulk of the work and is done in Section [HJ 
The properties of SRW needed are established in Section HI In Section we present 
the proofs of the main theorems. 

Our results here show some universality of the REM for dynamics of p-spin models 
with p > 3. This dynamic universality is close to the static universality of the 
REM, which shows that various features of the landscape of energies (that is of 
the Hamiltonian Hjq) are insensitive to correlations. This static universality in a 
microcanonical context has been introduced by [BM04J (see [BK06a, BK06bJ for 
rigorous results on spin-glasses). The static results closest to our dynamics question 
are given in [BGK06, BK07] where it is shown that the statistics of extreme values 
for the restriction of Hn to a random sets Xn C Sn are universal, for p > 3 and 
|-XjV"| = e c7V , for c small enough. 

This work was initiated during a concentration period on metastability and aging 
at the Max-Planck Institute for Mathematics in the Sciences in Leipzig. GBA and 
AB thank the MIP-MIS and Stefan Miiller for kind hospitality during this event. 
AB's research is supported in part by DFG in the Dutch-German Bilateral Research 
Group "Mathematics of Random Spatial Models from Physics and Biology" . 



2. Behavior the one-block sums 

In this section we analyze the distribution of the block-sums Y2i=i eie^^ v \ where 
Ci are mean-one i.i.d. exponential random variables, and {U{,% = 1, . . . , v} is a cen- 
tered Gaussian process with the covariance WJJJj = 1 — 2p\i — j\/N; v = v N is a 
function of N of the form 



v 



Liv-J. 



with uj E (1/2, 1). 



(2.1) 



As explained in the introduction, this process will serve as a local approximation of 
the corresponding block sums along a SRW trajectory. We characterize the distri- 
bution of the block-sums in terms of its Laplace transform 



E 



exp 



ue 



-■yN 



E 



ej.e 



(2.2) 
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Proposition 2.1. For all 7 such that 7 / '/3 2 G (0,1) there exists a constant, K = 
K(y, (3, u),p), such that, uniformly for u in compact subsets of[0, 00), 

lim tfi/V-^TW [1 _ r N {u)] = Kv?' p \ (2.3) 



Proof. We first compute the conditional expectation in (I2.2p given the a-algebra, U, 
generated by the Gaussian process U, 



E 



exp { - ue^ N Y e^^ } u] = 17 =r- 

i=i i=i 



^ i=i > 



(2.4) 



where 

g(x)=\n(l + x). (2.5) 

Note that importantly, g(x) is monotone increasing and non-negative for x G R + . 
We use the well-known fact (see e.g. |Sle61j ) that the random variables Ui can be 
expressed using a sequence of i.i.d. standard normal variables, Z i} as follows. Set 
Zi = {U x + U U )/{A - Ap{v - l)/Nf/ 2 and Z k = (U k - U k _ x ) /{Ap/N) 1 ' 2 , k = 2,...,u. 
Then Zi are i.i.d. standard normal and 

Ui = Y X Z X + ... + T l Z l - T i+1 Z i+1 - T U Z U , (2.6) 

where 

ri = ^1 -■£(!/- 1) and T 2 = ... = T U = ^. (2.7) 

Observe that Yli=i = 1- Let us define Gi(z) = Gi(z x , . . . , z v ) as 

Gi(z) = Vizi H h ViZi - T i+1 z i+ i Y v z v . (2.8) 

Using this notation we get 

1 = ^ (2^/2 ^ ? { — exp [ — ^ 3 (ue~^ N e 13 ^ *^ J I . (2.9) 

We divide the domain of integration into several parts according to which of the 
Gi(z) is maximal. Define D k = {z : G k (z) > Gi(z)Wi 7^ k}. On D k we use the 
substitution 



z i = b i + T i ( 1 N -logu)/(/3VN), ifi<h, 

' 2.10 
Zi = h - TifrN - \ogu)/(pVN), if % > k. 

It will be useful to define Yjj=i+i a j as ^2j=i a i ~ Sj=i a ji which is meaningful also 
for k < % + 1. Using this definition 

k 

G k (b)-G t (b) = 2 £ r„V (2.11) 

j=i+i 

Set 6> = — log(n)/(7iV) and define 

k j- 

D' k = {b: £ 6j + ^|*-i|(l + 0)>OVi^A;}. (2.12) 
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After a straightforward computation we find that (I2.9P equals 

-V 7 2 /2/3% (7 //3 2 f db -hYZ-,ti r -lVNG k (b)(l+9) 

(2.13) 



e " ' >" H u 

k=l 



X 



jl - exp ( - ^g ^e^^W-^^Wi^-^lfc-iKi+^j \ j 
i=i 



To finish the proof we have to show that u 1 ^ 2 is asymptotically the only dependence 
of (I2.13P on u (or on 9) and that the sum is of order uN~ 1 ^ 2 . We change variables 
once more to aj = bj/(l + 9) in order to remove the dependence of the integration 
domains on u. Then the sum (without the prefactor) in (12.131) can be expressed as 



£ 

k=l 



{i + eyda i (1+e)2Er=ig 2 

(2vr)<V 2 



(2-14) 

x |l - exp ( - ^g ^ e (^ G k( a )-^VP^U+i a ^- 2 P^ k - i \^ 1 + 9 ^ j | 

where D'l = {a : £* =i+1 % + ^f\k - i\ > Vi ^ fc}. 

Let 5 > be such that (1 + 5)7/7? 2 < 1, and let N > log(tt)/(7<5), so that \9\ < 6. 
We first examine the bracket in the above expression for a fixed k. On D" k 

V 

exp | _^g^ e (PVNG k (a)~2p^J2^ +1 a j -2p 1 \k-i\)(l+e)^ > gxp (.^(^(.JM)}. 
i=l 

(2.15) 

Write G k (a) as (recall (12TTD ) 

£ - cj log 

Crfc(a) = =. 2.16 



(2.17) 



The bracket of (I2.14p is then smaller than 

e -£( € - w]Qg ;v)(i + ^ 1 _ exp ( _ ^( e e-iogiV^| 

= ]v 2i ^ i e"^ (1+0) {l - exp ( - ^(eV^))}- 



The function e — exp ( — vg (e^/iA ) } is bounded for (el, uniformly in 

i/, if (1 + #h//3 2 < 1. Namely, if £ > 0, 



<f# (1+9) {1 - exp ( - vg (e*/u) )} < e~^ i+y) < 1. (2.18) 

If ^ < 0, then, since g(x) < x, 

{l-exp(-vg{e*/u))} < (l-exp(-e«)}, (2.19) 

which behaves like e^, as £ — > — oo. This compensates the exponentially growing 
prefactor, if (1 + 9)j / [3 2 < 1. Thus, under this condition, the bracket of (12.141) 
increases at most polynomially with N. 

In view of this at most polynomial increase, there exist 6 > small, such that the 
domain of integration in (12.141) may be restricted to Oj's satisfying 

V V 

v~ iy ^a 2 e (1 -6,1 + 6), h|<iV 1/4 , ^2\ai\ <v 1+s . (2.20) 

i=l i=l 
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The integral over the remaining aj's decays at least as e~ N& for some 5' > (by 
a simple large deviation argument). For all a satisfying (I2.20p . |Gfc(a)| < iV 1 / 4 + 
N~ 1 / 2 v 1+6 ' <C N 1 / 2 and thus, for any fixed u, uniformly in a, 

e -iVN GkW > l > and e -\TU^ >L (2 ' 21) 

Also, (1 + 0) v 1. Hence, up to a small error, we can remove all but the last 

occurrence of 6 in (I2.14p . 

Finally, taking X{ = ai for % > 2, xi = N 1/2 G k (a), and thus 

xi - Ap(x 2 H V x k - x k+1 x v ) 

a\ — , (Z.ZZ) 

(12.141) equals, up to a small error, 



dx e 2 ^'i=2 x i 



h 

x < I - 



x J D ,,T 1 Ny 2 (27ry/ 2 GXP I ~ ^ ~ 2ffivJ 6XP I " T + 2ffiV 
1 1 - exp - ^^^a+^ie-^^Wi^- 2 ^!^!)^)^ |. 



;2.23) 



i=l 



The last exponential term on the first line can be omitted. Indeed, 

-^ + ^ = Y ^[px 1 (x 2 + ----x u )-2 P 2 (x 2 + ..--x u ) 2 } ^0 (2.24) 

uniformly for all |xi| < N^ 1+s ^ 2 and \x 2 + ■ ■ ■ — x u \ < z/ 1+<5 )/ 2 , if S > sufficiently 

small. The integral over the remaining x is again at most e~ N . 
Now we estimate the integral over x 2 , . . . , x u , 

f exp ( - ^^eCi+^-e-^S^^^I^^a+^Y (2 . 25 ) 



where is the restriction of .D^' to the last v — \ coordinates (which does not depend 
on the value of the first one). Let V = (V 2 , . . . , V v ) be a sequence of i.i.d. standard 
normal random variables. Then, (12. 25ft equals 



VEDl 



V[V G S£]E[exp ( - J> ( e ( 1+ ^ I1 e-( 2 ^^Wi^+^|fc-i|)(i+^ 

i=l 

(2.26) 

The probability P[V G is bounded from below by the probability that the two- 
sided random walk, Ri = z^]=o Vj, i E 1*, with standard normal increments is larger 
than — jy/p\i\/fl for all i. This probability is positive and does not depend on N, 
which implies that, for all k, 

1 > P[V G D'l] > c> 0. (2.27) 

The expectation in (I2.26P is bounded by one, since the functions g is positive on 
the domain of integration. Moreover, as x\ — ► — oo, the argument of g in ( 12.261) tends 
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to zero (since the first exponential does, and the second is bounded by one on D'£) 
Hence 

(l+e)/^ -(2/VpE*U + i vsr+2p7|*-i|)(i+0)\ ^ M+e)pxi p-^v^E^U+i v j +2p rf \h-i\){i+ff) 



9 [e 
Therefore, as x 



-oo, 



(2.28) 



E 



cxp 



E 

i=i 



9 e 



p -(2/9VpEj= i+ i ^+2p 7 |fe-i|)(l+9) 



J _ e (l+e)/3a;i ]E r^ e -(2/3VpEj= l+ i^+2p7|fc-i|)(l- 



i=l 



3 (i+0)/3zi 



;i y^ E j" e -(2/3 v ^fl fc _ i +2 P7 |fc-i|)(i+e) 



Rk-i > 



7V^ 



(2.29) 



i=i 



/3 



— z| 



Since Ri is a centered normal random variable with variance a straightforward 
Gaussian calculation implies that 



E 



-(2(3 y /pR k ^+2p 1 \k-i\){l+e) 



Rk-i > 



k — i\ 



a 



- 7 2 p|fc-i|/(2/3 2 ) 



[2.30) 



Hence, (12.291) is essentially a summation of a geometrical sequence and therefore 
there exists constants c\, c 2 independent of k, such that 



1 - Cl e {1+e)(3xi < (J22SH < 1 - c 2 e (1 " 



■0)Pxi 



Vxi < 0. 



(2.31) 



Bounds (12.271) and (12.311) imply that (I2.25P is bounded from above and from below 
(with different constants) by 



CN~ 1 ' 2 



7 



6XP 1 " ^ " 2ffiV 



;i Aoe (i+0)/fe 1 ). 



(2.32) 



and hence (I2.23H is bounded from above and below by 



CuN- 1 ' 2 



jf danexp (- lx x - ^_)(1 A = M- 1 / 2 (2.33) 

Moreover, (12.251) is decreasing as function of min(fc, v — k). As this minimum tends 
to infinity, (I2.25P behaves as f{xi)N~ 1 ^ 2 which is of course satisfy the bound (I2.32p . 
Due to this convergence, the constants in the lower and the upper bound of (12.331) 
can be made arbitrarily close. This completes the proof of Proposition 12.11 □ 

We close this section with a short description of the shape of the valleys mentioned 
in the introduction. First, it follows from (I2.10p and the following computations that 
the most important contribution to the Laplace transform comes from realizations 
for which max{[/j : 1 < % < u} ~ 'y^/N / (3 with an error of order iV" 1 / 2 . It is the 
"geometrical" sequence in (I2.29P which shows that only finitely many neighbors of 
the maximum actually contribute to the Laplace transform. The same can be seen, 
at least heuristically, from a simple calculation 



E 



-c, 



P,7,P 



N 



;2.34) 
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Which means that, disregarding the fluctuations, the energy decreases linearly with 
the distance from the local maximum and thus the mean waiting times decrease 
exponentially. 



3. Comparison of the real and the block process 

We now come to the main task, the comparison of the clock-process sums with 
those in which the real Gaussian process is replaced by a simplified process. For a 
given realization, Y#, of the SRW, we set X%(i) = Hn(Yn(i)^ (the dependence on 
Yjv will be suppressed in the notation). Then X%(i) is a centered Gaussian process 
indexed by N with covariance matrix 

AJ = E[X° N (i)X N (j)} = R N (Y N (i) 1 Y N (j)) p . (3.1) 

Now we define the comparison process, Xjj(i), as the centered Gaussian process with 
the covariance matrix 

A^ n x^ N U }] J 1 -^-^ = IM. (3 ,) 

10, otherwise. 

For h G [0, 1] we define the interpolating process X^(i) = y/l — hX^(i) + \fhX\ii). 

Let £ G N, = to < ■ ■ ■ < U — T and U\, . . . , ut G M+ be fixed. For any Gaussian 
process X we define a function F N (X) = F N (X; {ti}, as 

e tkr(N) 

F N (X; {U}, { Ui }) = E[exp ( - £ £ ^*«) |*] (X) 

k=l i=t fe _ir(Af)+l 

3.3 

l t k r(N)-l _ 

= «p(-E E sO'* 1 " 

k=l i=t k _ ir (N) 

where r(N) = j\rV2 e iV7 2 /2/3 2 . Observe that E[F(X°;t,u)\y) is a joint Laplace trans- 
form of the distribution of the properly rescaled clock process at times ti. The 
following approximation is the crucial step of the proof. 

Proposition 3.1. If the assumptions of Theorem \l.l\ are satisfied, then for all se- 
quences {ti} and {ui}, 

jim^E^^^KK})^] -1^(4;^},^})] =0, y-a.s. (3.4) 

Proof. We use the well-known interpolation formula for functionals of two Gaussian 
processes due (probably) to Slepian and Kahane (see e.g. |LT91] 

tr(N) 



1 rl lr ^> r f) 2 F (X h ) l 

E[F N (X* N ) - F N (X° N )\y] = g I M E ( A & " A *H dX^dXti) \ y ] ■ (3 - 5) 

i¥=3 

We will show that the integral in (13.51) converges to 0. 
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Let k(i) be defined by tfc(i)_ir(iV) < i < tk(i) r {N). The second derivative in (13.51) 
is equal to 



=2-yN 



< 



U k (i)Uk(j)f3 2 N j3yr^ {X h {i)+X h {j)) 
e 2 7 N 



(3.6) 



x exp 



y * e 7_/V 



2^ 



where we used that g'(x) = (1 + x) _1 = exp(— g(x)) (recall (12. 5p ). and we omitted in 
the summation of F X (X X ) all terms different from i and j. To estimate the expected 
value of this expression we need the following technical lemma. 

Lemma 3.2. Let c G [—1, 1] and let U\, U 2 be two standard normal variables with 
the covariance K[UiUi] = c and A a small constant, < A < 1 — 7//3 2 (which will 
stay fixed). Define E N (c) = E N (c, P,j,u,v) and E N (c) = E N (c, fl,j,u,v, A) by 

uv/3 2 N i 



and 
E N (c 



-E 



exp \pVN{U 1 + U 2 )-2g( 



ue 



pVNUi-jN 



)-M 



ue 



f3VNU 2 -'yN 



) 

(3.7) 



(1-C)V2 eX P 



f _ 7 2 /V \ 
\ /3 2 (l+c) / ' 



l fc>( 1 /p 2 ) + \-l, 



/3 2 (l+c) J' ~ - \ HI- I ' " (3 ^ 

C'( 1 ,[3,u,v)Nexp{N((3 2 (l + c)-2 1 )}, if c < ( 7 //3 2 ) + A - 1, 

where C(j,/3, u,v, A) and C"(7, /3, u, v) are suitably chosen constants, independent of 
N and c. Then 

H/v(c) < E N (c). (3.9) 

Proof. Define k± = a/2(1 ± c). Let £7i, U 2 be two independent standard normal 
variables. Then U\ and U 2 can be written as 



(3.10) 



Hence, Ui + U 2 = k + Ui. Using g(x) + g(y) = g{x + y + xy) > g(x + y) and 
ue x + ve~ x > min(w, v )e' x ', we get 



»( 



ue 



f3VNU 1 ~'yN 



)+9( 



ve 



(3VNU 2 -'yN 



) 



> g I mm(u, v) exp ( h 



. 7 AT 

Denoting min(w,f) by u, we find that Sjy(c) is bounded from above by 
uv(3 2 N f dy 



(3.11) 



3 2 7 /V 



27 GXP 



2/1 + yl 



+ /3v / iV K+2/1 -2^( 



-K + Pv r Ny 1 /2+K-Pv r N\y 2 \/2- 1 N 



ue 



Substituting z\ — yi — PvNk + , z 2 = y 2 we get 
uv(3 2 N 2 K 2 +N/2 f dz_ 



o2yN 



27T 



exp 



z\ + z 2 



X 



exp ( — 2g\u exp < v^/V 



/3V 



(5k. 



>}■ 

(3.12) 



(3.13) 



\z 2 \ 
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The function exp(— 2g(ue^ x )) converges to the indicator function l x <o, as iV — ► oo. 
The role of x will be played by the bracket in the expression (13.131) . 

If this bracket remains negative for z close to zero, that is if 7 > —A' + (3 2 k 2 + /2 (or 
equivalently c < (•y/P 2 ) + A — 1), then the integral in (13.131) is bounded from above 
by 1. This yields the claim of the lemma for such c: 

3*(c) < ^^e^l N / 2 = C'( 1 ,P,u,v)Nexp (iV(/3 2 (l + c) -2 7 ) } = E N (c). (3.14) 



e 2 7 JV 

?2,,.2 



If this is not the case, that is 7 < —\'+/3k±/2, then we need another substitution, 

K-. . , T ( 27 



22 



f N 

V2 



vx -\v 2 \ - N[/3k. 



(3.15) 



'N 

This substitution transforms the domain where the bracket of ( 13 . 1 3f) is negative 
into the half-plain v\ < 0: The expression inside of the braces in (13. 13ft equals 
/3k + vi/2. Substituting ( 13. 15ft into (zf + zf)/2 produces an additional exponential 

prefactor exp ( — — -^s^r — )• Another prefactor iV -1 comes from the Jacobian. 

The remaining terms can be bounded from above by 

<lV exp t(/3 K+ -p-)(v 1 - -hi) - 2g(ue^ 2 )), (3.16) 



2n IV (3k + J \ k 

which can be separated into a product of two integrals. The integration over v 2 gives 
a factor 

((/?«+ - ^)^rj l < C{\)k-J < C(A)(1 - c)- 1 ' 2 . (3.17) 

Using properties of g, the integrand of (13. 16f) behaves as exp{— 2v ij//3k + } as v\ — > 
00, and as exp{(/5ft + — (2 7 //3k + ))ui} as v% — > —00. Therefore, the integral over v\ is 
bounded uniformly by some A-dependent constant for all values of c > — 1 + (7 / ' f3 2 ) + 
A. Putting everything together 

*„(c) < C(l - cl-^^e^^iexp ( - ^± = 2l?N 
~N{c) cj ^ e ^exp^ 



2/3 2 /^ 

L7( 7 ,/5,M,t;,A)(l -c)- 1/2 exp<| - - J, - f = -\ (<■). 



(3.18) 

■ X = ~, r (r\ 

(3 2 {l + c) 

This finishes the proof of Lemma 13.21 □ 

Let \\d\\ = min(d, N — d) and = dist(Y N (i) ,Y N (j)) . Define, with a slight abuse 
of notation, A° = (1 — 2dN~ l ) p . That is K° d is the covariance of X^(i) and X^{j) 
if Dij = d. The next proposition, which will be proved in Section HI will be used to 

L 7V" 



control the correlations of the process X° 



Proposition 3.3. Let 7 and /3 satisfy the hypothesis of Theorem ! 1. 11 and let v be as 

in (12.11) . Taen, /or any 77 > 0, taere exists a constant, C = C(/3,j, u, rj), such that, 
y-a.s. for N large enough, for all d 6 {0, ... , N} 

tr(N) r /aa 1 

^ = d} < C tV^^r" j +tr{N)v- 1 e^ , (3.19) 
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tr(N) 



E i{^ = rf}(AS-Ai)<^^l{d<i/}. 



i,j=l,ij^j 
[i/v]=\j/v] 



(3.20) 



We now conclude the proof of Proposition 13. 11 that is we prove that the right-hand 
side of (13.51) tends to 0. Observe first that Dij is smaller than \i — j\. Hence, for 

A° ij =[l-2N- 1 D ij ] p >[l-2N- 1 \i-j\Y>A 1 ij . (3.21) 

Since Ajj = for with [i/u\ ^ A?- - A-,- < if and only if A°- < 0. 

The summands on the right-hand side of (13.51) can be written as differences of two 
non-negative terms: 



(A 



Aj,-)+E 



d 2 F N (X 



NJ 



y 



(A°)-E 



d 2 F N (X 



AO 



dX(i)dX(j) 



y 



(3.22) 



.dX(i)dX(j) 

We bound this expression using Lemma I3.2[ For given {ui} let 

Sat(c) = max{S A r(c,/?,7,M i ,M i ) :l<i,j< £}. (3.23) 

Then Hjv(c) satisfies (13.81) for some constants C and C" and it is therefore increasing 
in c. The absolute value of the right-hand side of (13. 5p is then bounded from above 
by 

MAO _ o*} i — i / t/"D \ . HN) 



E(4-4)+ E 



&F N {X° 



NJ 



dX(i)dX(j) 



N ( tr(N) -, 

^E E 1{A, = 4(A2) + / 



^ +E(4)- E 



<9 2 i^(Al 



&3 



dX{i)dX{j)\ 



E(hA»)dh 



\i/Zft\j/»} 



tr(N) 

E HD 



(3.24) 



d}(AS-Ai)H(Ag 



\i/"] = Wv] 



tr(N) 

+ E l{^- = d}(AS)-S(0) 

»j':|»-j|>JV/2 

From the definition of H it follows that, 



jf' =(Md A < Cexp { - ^g-j} jT'd - M-v 2d ,, 



(3.25) 



The last integral can be easily evaluated and is smaller than 2 for all c G [—1, 1]. 
Using Proposition 13.31 the first line of (I3.24p is smaller than the sum of the following 
two terms: 



cf>(A0 2 2^Wp{ 



j 2 N 



and 



N 



tr(N)e 



n\\d\\ 



d=0 



V 



iP(l + A») 



(3.26) 



(3.27) 
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The second line of (13.241) is bounded by 

tr(N)d 2 ~ 4flN 

d=0 

The third line is non-zero only if p is odd, and in that case it is bounded by 

J2 C \t\(N) 2 2- N f N d ) + tr(N)u-^] - l) P H(0), (3.29) 

d=N/2 L ^ ' -I 

We estimate (I3.26P first. Let I{u) be defined by 

I(u) = u log u + (1 - u) log(l - u) + log 2, (3.30) 

and let 

Stirling's formula yields Jn(u) > (4w(l — w)) _1 uniformly in u on compact subsets 

of (0, 1). Further, Jn(u) < CN 1 ^ 2 for all u G [0, 1]. From the definitions of r(N) and 
S, we find that 

mB = cjr t 2 N^ 2 (l - 2 ^) P exp { JVT pA7 (|) } Jv (|) , (3.32) 

where 

T (u) = tf-W- nug-W if (1-2«)p>^ + A-1, 

P,/3M J \ J - /(«) + /3 2 (1 + (1 - 2m)p ) - 2 7 , if (1 - 2uf < $ + A - 1. 

(3.33) 

Lemma 3.4. There exists a function suc/j that for allp > 2, and 7, /? satisfying 
7 < C(p)P an d 1 < P 2 > there exist positive constants 5, 5' and c such that 

T pAy (u)<-8 /or all u G [0,1] \ (1/2-5', 1/2 + 5'), (3.34) 

and 

T Pl/ 3, 7 (w) < -c(u - 1/2) 2 for all ue (1/2 -6', 1/2 + 6*). (3.35) 
Moreover ((p) is increasing and satisfies (11.101) . i/iai 

C(2) = 2~ 1/2 , C(3) = 1.0291, and lim CO) = V 2 lo g 2 - ( 3 - 36 ) 

Proof. Since 7//? 2 < 1, the second line of the definition of T Pi/ 3 i7 is used only for 
p odd and u > u c (p, @, j, A) = (1 + (1 - A - 7//3 2 ) 1 / p )/2 > 1/2. Furthermore, 
T pA7 (l/2) = V pA7 (l/2) = and 

T" (1/2) = J 4 (^" ~ X )' lfp = 2 ' (3.37) 
p ' ' 7 1 —4 otherwise. 

The second derivative is always negative for /3, 7, p satisfying the assumptions of 
Theorem ll.il Therefore (I3.35f) holds. 

The second line of the definition of T P)( g )7 (u) is decreasing in u. Hence for u > u c 

T pA7 (n) < T pA7 ( Mc ) = - 7 (1 - 1 /(3 2 ) - I(u c ) (3.38) 

which is obviously strictly negative and H3.34I) is proved for u > u c . 
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Figure 1. Function T Pi7;/ g for p = 2, 3, 4 and various values of j/P- 

For any 5' > and u < 1/2 — 5' the function I(u) is strictly positive, and the 
function = 1 — 1/(1 + (1 — 2u) p ) is bounded. Therefore, if 7//? is sufficiently 
small, then T p ^ a {u) < —5. If p is even, the function T P)j a i7 is symmetric around 
u = 1/2. If 1/2 < u < u c (p, (3, 7) and p is odd, then 

T pA » < T Pilj0 (n) = -I(u) < (3.39) 

and the proof of (I3.35P is finished. 

To prove the first part of (13.361) we should check that (13.351) holds for all 7 < 
2 _1 / 2 /3. However, T 2j g 7 (u) is increasing in 7 2 //3 2 and I[u) > (1 — 2u) 2 /2. Thus, for 

1 < 2- 1/2 /3, 

^-M<-\( l - l + { l-2u T ) -\ {1 - 2 ^ < X40 > 
The right-hand side of the last inequality is equal for u = 1/2 and its derivative 

2(l-2n)(l- (1 + (1 1 _ 2M)2)2 )>0 for all u< 1/2. (3.41) 

The symmetry of T2„a i7 around 1/2 then implies the first part of ( 13.36!) . 

Obviously, $(0) = 1/2, $'(0) = -2p, 1(0) = log2 and J'(0) = -00. Hence, 
for j/p = A/log 2 there exists w small such that T Pt p n (u) is positive. This implies 
C(p) < V2 log 2. If w e (0, 1/2) then limp^ $(m) = 0. This yield the second half of 
(E3E1). 

For illustration you find the graphs of function T Pi ^ j7 for p = 2,3,4, (3 = 1, 
and 7 = (solid lines), 7 = \flj2 (dashed lines), 7 = 1 (dash-dotted lines) and 
7 = ^/2 log 2 (dotted lines) on Figure [TJ The value of ((3) was calculated numerically 
using the figure for p = 3. □ 

We can now finish the bound on ( I3.26p . Lemma [3.41 and bounds on the function 
Jn yield that for d/N ^ (1/2 — 5', 1/2 + 5') the summands decrease exponentially in 
N. Therefore they can be neglected. The remaining part can be bounded by 

(l/2+«')iV 

C Yl t2ArV2 ( 1 -|) exp(-ciV( C //iV-l/2) 2 ) 



d=(l/2-6')N 

-6' 

< Ct 2 N z ' 2 I x p e~ c ' Nx2 dx 



(3.42) 



5' 

00 

,/„2 . N^oo 



if p > 3. 
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Similarly, for (13.271) we have 

N/2 



(EOTI) < CV# 2 i/ 4 (l - -^j exp{NT{d/N)), (3.43) 



2d\p 

tl\*'"i> '11- 

d=0 

where, setting = min(u, 1 — u 



j , > J 2/3 2 /3 2 (1+(1-2w)p) + ^H^ll' if ( 1 2m) p > J 2 + A 1, 

^ + /5 2 (1 + (1-2 M )p)-2 7 +77|H|, if (1-2u)p< $ + 



(3.44) 

Observe first that the second part of the definition of T Pi ^ i7 is always strictly negative. 
It is also easy to be checked that it is possible to choose 6, 5' and r\ small such that 
the first part of the definition of T(w) < S for all ||w|| > 5'. Therefore such d can be 
neglected. Around d = the function T(x) can be approximated by a linear function 
— cx, c > 0, and the summation by an integration. As an upper bound we get 

CtN 3 ' 2 v~ l f e~ cNx dx < CtN 1 / 2 ^ 1 0. (3.45) 

Jo 

An analogous bound works for d close to N and p even. 
For (I3~28|) we have 

V 

(E2ED < C^tN~ 3/2 d 2 [l - (1 - 2dN- 1 ) p }~ 1/2 exp{NT{d/N)). (3.46) 

d=0 

The linear approximation of T and of the bracket in the last expression yields an 
upper bound 

CtN 3 / 2 f x 3 ^ 2 e- c ' Nx dx < CtN- 1 ^=^> 0. (3.47) 
Jo 

Finally, since S(0) = Ce~ Nl2 ^ 2 , it is easy to see that the second half of (13.291) 
tends to 0. The first half equals (up to constant) 



d=N/2 



<Ct>{ £ *W^ + £(^-l)V/V*/»} 



d>N/2+N 3 / 5 



where we used the known approximation of (^) < CiV 1 / 2 2 N e * 2 / 2Ar for d = (N+i)/2 
and z <C iV 2 / 3 . The first term in ( 13.481) tends to by a standard moderate deviation 
argument. The second one can be approximated by 

POO 

Ct 2 N i-(j>/2) / x P e -* 2 /2 dx o (3.49) 

Jo 

for p > 3. This completes the proof of Proposition 13.11 □ 
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4. Random walk properties 



In this section we prove Proposition 13.31 For A C Sn let Ta = min{&; > 1 : 
Y/v(fc) G ^4} be the hitting time of A. We write ¥ x for the law of the simple random 
walk Y N conditioned on Y N (0) = x. Let Q = Qi, i G N, be a birth-death process on 
{0, . . . , iV} with transition probabilities Pi 7 i-i = 1— Pi,i+i = i/N. We use Pk and Ek to 
denote the law of (the expectation with respect to) Q conditioned on Q Q = k. Under 
Poj Qi has the same law as dist(Y/v(0), Yn(i)). Define Tk = min{z > 1 : Qi = k} the 
hitting time of k by Q. It is well-known fact that for k < I < m 

Pi[T m <T k ] = irr_i J ,'-r (4- 1 ) 



Em-i /n-i\ 
i=k 



(V) 



Finally, let Pk(d) = Po{Qk = d). We need the following lemma for estimating pk(d) 
for large k. 

Lemma 4.1. There exists K large enough such that for all k > KN 2 log iV =: JC(N) 
and x, y 6 Sjy 

F y [Y N (k) =xUY N (k + l) 



x\ 



and thus 



Pk{d) + Pk+i{d) _ N f N 
2 [d 



-N 



< 2 



< 2 



-AN 



-8N 



(4.2) 



(4.3) 



Proof. The beginning of the argument is the same as in [Mat87] . We construct 
coupling between Y/v (which by definition starts at site 1 = (1, . . . , 1) e Sn) and 
another process Yjy. This process is a simple random walk on Sn with the initial 
distribution [i* N being uniform on those x G Sn with dist(x, 1) even. The coupling 
is the same as in |Mat87j . This coupling gives certain random time T~n which can be 
used to bound the variational distance between /i* and the distribution p^ of Y N (k): 
for k even 

tU/4,/4) = max \n* N (A) - p k N {A)\ < F[T N > k}. (4.4) 

The law of Tn is as follows. Let U = dist(Y^(0), 1). That is U is a binomial random 
variable with parameters iV and 1/2 conditioned on being even. Consider another 
simple random walk Yjj on Su started from 1. The distribution of is then the 
same as the distribution of the hitting time of {x G Su : dist(l,x) = U/2}. It is 
proved in [Mat87j that P(T N > NhgN) — ► c < 1. It is then easy to see that, 

W[Tn > £(N)} < c KN/2 < 2~ m } (4.5) 

if K is large enough. Thus, for even k > JC(N), docfax, p k N ) < 2~ 8N and thus 
\p%(x) — p k N (x)\ < 2~ 8N for all x G <Sjv- A similar claim for k odd is then not difficult 
to prove. The second part of the lemma is a direct consequence of the first part. □ 

Lemma 4.2. Let 7, (3, v satisfy the hypothesis of Proposition lSUl Then, there exists 
a constant, C = C(/3, 7, v), such that for all N large enough, y-a.s. 

tr(N) 

!{Aj = d}< Ctr(N)l{d < u}, (4.6) 

[i/u\ = [j/u\ 

and for all d G {0, ... , N}. 
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Proof. The lemma is trivially true for d > v. For d < u, let 

V 

p(d) = E J2HQi = d}. 



(4.7) 



i=l 



We have p(0) > N 1 and p(d) > Po[T d < v\. This probability is decreasing in d and 

N N - 1 N-u+1 



> e 



-v 2 /N 



N N N 
Thus p(d) > e~ v2 l N for all d < v. To get an upper bound on p{d) we write 



p(d) <E [J2 HQi = d}]=l + E d [j2 1 iQi = d} 

i=l i=l 



1 + 



1 



Pd[T v < T d ] 



(4.8) 



(4.9) 



However, using (14.11) . 

N — d N — d ( N ^\ 1 
P d [T v < T d ] = -^P d+1 [T V < T d ] = - t = 1 - 0(vN-'). (4.10) 



N 



N 



Eu-l (N—l\ 
i=d { i ) 



Since v < N, p(d) < 2. 

Consider now one-block contribution to (14.61) 



J2 l{D {j = d} =: v 2 Z. 



(4.11) 



Of course, Z e [0, 1] and, using the results of the previous paragraph, 

e- v2/N (2v)- 2 < E[Z] < 2v~\ (4.12) 

The left-hand side of ( 14.61) is stochastically smaller than v 2 Y^k=i where Zk are 
i.i.d. copies of Z and m = \tr{N)/v\. By Hoeffding's inequality [Hoc63j, 

m 

P[^Z fe > 2mE[Z k ] < exp{-2m 2 E[Z k } 2 } < exp{-mV 2iy2/7V (2z/T 4 }, (4.13) 
i=i 

where we used the lower bound from (14. 12ft . Since u/N 2 <C N, by the Borel-Cantelli 
lemma, the left-hand side of ( 14. 6 ft is a.s. bounded by 

u 2 2mE[Z} < Ctr(N) (4.14) 

for all N large enough and d < v. This completes the proof of Lemma 14.21 □ 

Proof of Proposition ^. 3\ We prove (I3.20f) first. Observe that for i,j in the same 
block 

2p\i-j\\ 2p(\i-j\-d) . „{d 2 \ 



A - A 1 



2d\P 
1 ~ iV\ 



N 



N 



+ [w ). (4.15) 



The contribution of the error term is smaller than the right-hand side of (I3.20p . as 
follows from Lemma 14.21 

To compute the contribution of the main term, let 



V 

p(d) = E [J2(i-d)l{Qi = d} 

i=i 



(4.16) 
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Let T\ = T d and T| = min{i > T%~ 1 : Q { = d}. Then 

oo oo 

p(d) = E [J2(n - d)l{Tl < u}] = E [ - Tl + Tj - d)l{T> < u} 



3=1 3=1 

oo 

< E [(T d - d)i{T d < u}\ (i + unun <v-d}\ 



i=l 



(4.17) 



Using gSD, P [T d = d]< Ce~ s l N and further 



Po[Ti > d + 2 k] <(^ k )(±y<Cg. (4,8) 



Hence, cd 2 N~ l < E [(T d - d)l{T d < u}] < Cd 2 N~\ 
For the second term in (14.171) we write 



l+J2E d [T^l{Tj<u-d}) 

oo 

< 1 + E d [T d l{T d < v - d}} (l + E d [T l d l{r d <u- d}]) (4.19) 

i=i 

oo 

= J2{E d [TdHT d <u-d}]} k . 



Using the well-known estimate ( 2 fc fc ) < Ck 1 / 2 2 k and k <2 k , 

E d [T d l{T d < f>(f) (£)* < C± (£)' < (4.20) 

k=l V / k=l 

and glSD is finite. Thus < Cd 2 N~ 1 for all d e {0, . . . , v). 

The one-block contribution of the first term of (I4.15P to (13.201) is then given by 

| I>'-il - d) 1 ^ = rf > = : f ^ ( 4 - 21 ) 

«,i=i 

with Z e [0, 1] and 

cd 2 N- l v^ < E[Z] < Cd 2 N- l u- 2 . (4.22) 

Therefore, as in the proof of Lemma l4~2l Hoeffding's inequality and (14.221) imply that 
the contribution of the first term of (14. 1 5j) to (13.201) is smaller than Ctr(N)d 2 N~ 2 , 
which was to be shown. 

Finally, we prove H3. 19D . Since we are interested in an upper bound only we can, 
without loss of generality, restrict the summation on i < j. We first consider the 
contribution of pairs (i, j) such that j —i > K,(N). Then necessarily, \i/n\ ^ \J/n\. 
Let R = tr(n). Lemma [4 . 1 1 yields 



E[ 1 iD ij = d}]= P3-i(d)<CR 2 2- N (^\ (4.23) 

j-i>K{N) 3-i>fC(N) V ' ' 
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Further, 

r R 

Var HDij = 4 

'-j-i>K(N) 

= E E P [Au* = D *** = d ~\ -*l D h* = d M D n,h = d] ■ 

h-h>K(N) j 2 -i 2 >K.(N) 

We can again suppose that %\ < i 2 . The right-hand side of (14.241) is non-null only 
if i\ < %2 < ji < j 2 or i\ < %i < ]2 < ji- We will consider only the first case. The 
second one can be treated analogously. In is not difficult to see using Lemma 14.11 
that if %2 — ij > /C(iV) or j 2 — ji > then the difference of probabilities in the 

above summation is at most 2~ 4N . Therefore, the contribution of such (ii, i 2 , ji, j 2 ) 
to the variance is at most i? 4 2~ 4Ar . 

If %2 — i\ < K.{N) and j 2 — ji < then, using Lemma fl~T1 again. 

F[D ildl = A 2J2 = d]< C2- N f^j . (4.25) 

We choose e > 0. For \\d\\ < (1 — e)N/2 we have 



E E V[D iltjl = D i2d2 = d] 



h-h>K{N) j 2 -i 2 >K(N) 
i 2 -h<K{N) j 2 -h<K(N) 



< CK{N) 2 R 2 2~ N ^ < CJC(N) 2 R 2 e- NI ^- e ^^ < N^R 2 u~\ 

(4.26) 

say. For > (1 — e)N/2, that is \d — N/2\ < eN/2, we have for e small enough 
(how small depend on 7 and (3) that 2~ N '(*[) > N 7 R~ 2 . Then, 



E E V[Dn, n =D l2 , 2 = d] 



ji-ii>K(N) j 2 -i 2 >tC(N) 
h-ii<K{N) j 2 -j x <K{N) 



(4.27) 



We have thus found that the expectation of the summation over j — i > JC(N) is 
smaller than the right-hand side of (I3.19P and the variance of the same summation 
is much smaller than N~ 3 times the right-hand side of (13.191) squared. A straight- 
forward application of the Chebyshev inequality and the Borel-Cantelli Lemma then 
gives the desired a.s. bound for pairs j — i > JC(N) and all d 6 {0, ... , iV}. 

Choose again e > 0. For j — i < IC(N), observe first that if ||<i|| > (log N) 1+e 3> 
logiV then the summation over such pairs in ( 13 . 1 9f) is always smaller than 

)C(N)R < Rv- l e^ for all 77 > 0. For the remaining d% that is < (logiV) 1+e ', 
let Kn > K be the smallest constant such that K^N 2 log N is a multiple of v. Since 
v <C N 2 , K N — K <C 1. As the difference between K and K N is negligible, we will 
use the same notation JC(N) for K^N 2 log iV and we will simply suppose that /C(iV) 
is a multiple of v. The summation in (13.191) for j — i < K.(N) can be bounded from 
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above by 

tr(N) K(N)-1 \R/K(N)~\ K{N) 

,K{N)i+k+m — d}, (4.28) 

Q<j-i<K(N) k=0 E=0 m=j k 

where jk is the smallest integer such that [(K.(N)£ + k)/u\ ^ \ {K,(N)i + k + jk)/v\, 
which does not depend on I. We define random variables Zt(J, d) by 

K{N) 

z e(j, d) = y[m) E 1 { D ic(N)e+k,K(N)i+k+m = d}. (4.29) 

^ ' m=j 

The sequence {Ze(j, d) : i > 0} for fixed j and g? is a sequence of i.i.d. variables with 
values in [0,1]. 

Let E N = {d: \\d\\ < (\ogN) l+e \d> N/2}. For d e E N 

,(M) > 0] < (^)^(Ti < fc(JV)) < Qe^[e-^/^ 10 ^]. (4.30) 



According to Lemma 3.4 of [CG06] . 

E d [exp(-AT 1 m(iV)- 1 )] < (2~ N m(N)X~ 1 + £ N (d))(l + o(l)), (4.31) 

for iVlogiV « m(iV) « 2", with £ n {h) = 2""f Q~ a ( fe ^.)i. Taking m(iV) = 

N 2 and d G .Eat it is not difficult to check that for e small enough 

E Zd [e~ XT ^ N2 ] <2- N( - 1 ~ £ l (4.32) 

Hence, 

K{N)-1 \R/K(N)] 

'[ U { E E Ze(j k ,d)>0}] 



P 



/ ]\T \ 

i-e'N 



- c (r(io g ^)-i) fl(logiV)1+ ' 2 " < " > - c ' r ''' 

(4.33) 

for some e' small. Hence, d G En do not pose any problem, by the Borel-Cantelli 
lemma again. 

To treat d < (logiV) 1+£ we will distinguish two cases: jk < 2d and jk > 2d. 
For the first case, observe that for any d < v there are at most dK,{N)/v values of 
k G {0, . . . , K.(N) — 1} such that j k < d. Further, as before, Zi(j k ,d) < Zi(0,d), 
E[Z e (0,d)] > 1/(JV/C(JV)), and E[Z t {0,d)] < C/K(N). Hence, by Hoeffding's in- 
equality, the probability 



P 



\R/K[NJ\ 



decreases at least exponentially with N and thus for jk < 2d, a.s, 

jc(n) z ^ d ^Ym- ( 4 - 35 ) 
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For j > 2d and N large enough, Zg(j, d) < Z e (d + 6, d). We have, 

cN~ 6 < K(N)E[Z e (d + 6, d)} < CN~ 3 . (4.36) 

Indeed, the lower bound is trivial and for the upper bound we use the fact that the 
probability that Yjv reaches d before returning to d + 6 is smaller than CN~ 5 and 
before the time 1C(N) there are at most JC(N) tries. Hence, for j > 2d the probability 

P[/C(A0 £ Z e (k,d) > j^^} (4.37) 

decreases at least exponentially in N and thus the interior inequality is not valid 
a.s. for all N large. Summing over k we get 

K(N)-1 \R/!C(N)] 

(4.38) 

since 7//5 2 < 1. □ 



5. Convergence of clock process 



We will prove the convergence of the rescaled clock process to the stable sub- 
ordinator on space D([0, T],M) equipped with the Skorokhod Mx-topology. This 
topology is not commonly used in the literature, therefore we shortly recall some of 
its properties and compare it with the more standard Skorokhod Ji-topology, which 
we will need later, too. For more details the reader is referred to [Whi02j for both 
topologies and to [Bil68j for detailed account on Ji-topology. 

5.1. Topologies on the Skorokhod space. Consider space D = D([0, T],M) of 
cadlag functions. The J\ -topology is the topology given by the Ji-metric: for f,geD 

djAf,9) = } ni A\\f oA-fi-IU V || A — e||oo}, (5.1) 
AeA 

where A is the set of strictly increasing functions mapping [0, T] onto itself such that 
both A and its inverse are continuous, and e is the identity map on [0, T]. 

Also the Mx-topology is given by a metric. For / G D let Tf be its completed 
graph, 

r J = {MGlx[0,T]:z = «/(*-) + (1 - a)f(t), a G [0, 1]}. (5.2) 

A parametric representation of the completed graph Tf (or of /) is a continuous 
bijective mapping <j>(s) = (0i(s), (j) 2 (s)), [0,1] t— > Tf whose first coordinate is 
increasing. If n(/) is set of all parametric representation of /, then the Mi-metric 
is defined by 

d Ml (f,g) = MiWfa - ih\\oo V - ^2||oo : <P e n(/), tp G n(^)}. (5.3) 

The space D equipped with both M\- and Ji-topologies is Polish. The Mx-topology 
is weaker than the Ji-topology: As an example, consider the sequence 

/„ = 1{[1- l/n,l)} + 2 -1{[1,T]}, (5.4) 

which converges to / = 2 • 1{[1,T]} in the Mx-topology but not in the Jx-topology. 
One often says that the Mx-topology allows "intermediate jumps". 
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We will need a criterion for tightness of probability measures onD. To this end 
we define several moduli of continuity, 

w f {6) = sup { min - /(Ol, \f(t 2 ) - f{t)\) : h < t < t 2 < T,t 2 - t x < 8}, 

w'(8) = sup{ inf If^-iafitJ + il-^fih))] : h < t < t 2 < T,t 2 - t x < 8}, 

v f (t, 8) = sup {!/(*!) - f(t 2 ) I : h, t 2 e [0, T] U (t - S,t + 8)}. 



(5.5) 



The following result is a restatement of Theorem 12.12.3 of [Whi02] and Theorem 15.3 
of pl68] . 

Theorem 5.1. The sequence of probability measures {P n } on _D([0,T],M) is tight in 
the Ji-topology if 

(i) For each positive e there exist c such that 



Pn[f: 



> c < e, 



n > 1. 



(5.6) 



(ii) For each e > and 77 > ; there exist a 8, < 8 < T, and an integer Uq such 
that 



Pn[f : W f (8) > V ]<S, 



n > n , 



(5.7) 



and 



P n [f :v f (0,8) >rj\ <e and P n [f : v f (T,8) > 77] < e, n > n . (51 
The same claim hold for the Mi-topology with Wf(8) in (15.71) replaced by w'f(8). 



5.2. Proof of Theorem ll.il To prove the convergence of the rescaled clock process 
Sn{-) = e~ lN S]y(-r(N)) to the stable subordinator V^/pi, we check first the conver- 
gence of finite-dimensional marginals. As can be guessed, Proposition 13.11 will serve 
to this purpose. Let £, {ui} and be as above. Then, 



E 



exp I -^2ui(S N {t k ) - 5^(4-1)) 



i=l 



Yi 



N 



(5.9) 



= E[F N (X° N ;{t l },{u l })\Y N ] = E[F N (Xlr;{U}, {m})] + o(l), 

as follows from Proposition 13.11 

The value of E[F/v(X^; {U}, {««})] is not difficult to calculate. Define jn{ 
[tir(N)/v\. Then 



E[F 7V (X^;{t i },{ Mi })] =E 



> E 



[n n «p 

k=lj=j(k-l)+l 



exp 



Uk 



t k r(N)-l 
k=l i=t k _ ir (N) 



u-l 
i=0 



(5.10) 



Since the process is a piece-wise independent process, the product in ( 15.101) is a 
product of independent random variables. Then expectations of all of them can be 
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then bounded using Proposition 12. 11 We get, for 5 > fixed and N large enough, 

i jjv(fc) 

E[i^p& ; {uh k»] > n n 

fe=ii=w(fe-i)+i 

>JJ (!-(! + ^z/JV-^e-^/^^/^j^W-i^C*-!)-! (5.11) 
fe=i 

> JJexp { - (1 + 2(J)(t fc - t fc _!)if^}, 
fc=i 

which is (up to 1 + 25 term) the Laplace transform of V^/^^K-). A corresponding 
upper bound can be constructed analogously. 

To check the tightness for S N in D([0, T],R) equipped with the Skorokhod Mi- 
topology we use Theorem 15. 11 Since the processes Sn are increasing, it is easy to see 
that condition (i) is equivalent to the tightness of the distribution of Sn{T), which 
can be checked easily from the convergence of the Laplace transform of the marginal 
at time T (the limiting Laplace transform tends to 1 as u — > 0). 

In order to check condition (ii), remark that for increasing functions the oscillation 
function w'g (8) is always equal to zero. So checking (ii) boils down to controlling 
the boundary oscillations ^,§^(0, 5) and v§ (T, 5). For the first quantity (using again 
the monotonicity of Sn) this amounts to check that P[Sat(<5) > rf\ < e if 5 is small 
enough and iV large enough. Using the convergence of of marginal at time 5, it is 
sufficient to take 5 such that F[Vy/p2(K5) > rf\ < e/2, and take n such that for all 
n > no 

\nS N (6) >v\- nV 7/f3 2(K5) >T)}\< e/2. (5.12) 
The reasoning for v§ N (T, 5) is analogous. □ 

5.3. Coarse-grained clock process. To prove our aging result, that is Theo- 
rem II. 2\ we need to modify the result of Theorem 11.11 slightly. Let Sn be the 
"coarse-grained" clock processes, 

S N (t) = ^ w S N (u[tr(N)u- 1 \). (5.13) 

For these processes we can strengthen the topology used in Theorem that is we 
can replace the Mi- by the Ji-topology. 

Theorem 5.2. If the hypothesis of Theorem \l.l\ is satisfied, then 

S N (t) V l//3 2(Kt) y - a.s., (5.14) 

weakly in the Ji-topology on the space of cadlag functions D([0,T],R). 

Unfortunately, we cannot prove the theorem with estimates we have already at 
disposition. We should return back and improve some of them. First we show that 
traps with energies "much smaller" than 'y^/N / (3 almost do not contribute to the 
clock process. Let B m = 7-y/iV / (3 — m/ ((3\/N) and let 

Ltr(JV)J 

S%(t) = e^ N e * ex P < B m }- (5.15) 

i=0 
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Lemma 5.3. For every T and rj, e > there exists m large enough such that 

P[S£(T) > V \y] < e, y-a.s. (5.16) 

Proof. To prove this lemma we should improve /modify slightly the calculations of 
Sections [5] and El With the notation of Section [5] define 



exp 



{ " ^ N E ^ Ui m < B m }}} . (5.17) 

1=1 



(comparing with (12.21) observe that we set u — 1). We will show that 

lim f(N)e N ^ 2 ^ 2 [1 - = K m , {5. It 



N— >oo 



with K m — > as m — > oo. The proof of this claim is completely analogous to the 
proof of Proposition 12.11 One should only modify the domains of integrations. More 
precisely, the definition of D k which appears after (12.91) should be replaced by D™ = 
D k n{z : G k (z) < B m }. Hence, D' k becomes D' k m = D' k n{b : G k {b) < -m/(J3/y/N)}, 
which then restricts the domain of integration in ( 12. 33ft to (— oo, —m//3]. Hence, the 
constant K m can be made arbitrarily small by choosing m large. 
Further, as in Section [31 define 



Tr(N)-l 

Y '(X) = exp(- E 9(e- rN e PVWX(i) HX(i)<B m })). (5.19) 

i=0 



Then, as in Proposition 13.11 we will show 

lim E[F%(X° N )\y] —E[F%(Xlf)] = 0, >a.s. (5.20) 



We use again (13. 5p to show this claim. Although the indicator function is not dif- 
ferentiable, we will proceed as if it was, setting (l{x < B})' = —5(x — M), where 5 
denotes the Dirac delta function. As usual, this can be justified e.g. by using smooth 
approximations of the indicator function. The second derivative of F™(X) equals 



■ 2 f_ 

e 2yN 
X 



(3VN / V pVN 
< u *p* Ne PVN(xkM+xW))-*rN exp ( _ 2g (ue^ xh ^ N ) - 2g(ue^ x ^^ N )) 

x < B m} - S -^l) (l { X U) < B m} - ^Wll) . 

(5.21; 



We should now bound the contributions of four terms. The one with the product of 
two indicator functions is easy, because we can use directly the result of Lemma 13.21 
For remaining three terms, those with the product of one indicator and one delta 
function, and this with two delta function, the calculation should be repeated. How- 
ever, in the end we find that (I5.2ip is bounded by E(Gov(X(i),X(j))) as before. The 
presence of the delta functions makes actually the calculations slightly less compli- 
cated. The proof then proceed as in Section [3l 
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We can now finish the proof of Lemma 15 .31 By (15.171) and (15.201) , 
E[exp(-^(T))|y] =E[F%(X° N )\y] =E[F%{Xj r )\y] +o(l) 

= (1 - # m /(iV)-V JV T' 2 / 2/32 ) rr(A0A/ + o(l) = e~ KmT + o(l)- 
Since K m — > as m — > oo, 



(5.22) 



mm > , M < LMSH (5 . 23) 

can be made arbitrarily small by taking m large enough. □ 

We study now how the blocks where the process visits sites with energies larger 
than B m are distributed along the trajectory. To this end we set for any Gaussian 
process X 

s%(i;X) = l{3j :iu<j<(i + l)u,X(j) > B m }. (5.24) 
and we define point process H^{X) on [0,T] by 

Tr{N)/u 

H%(X;dx) = s™(i;X)5 il//r{N) (dx). (5.25) 

Lemma 5.4. For every m G M the point processes H^(X Q N ) converge to a homoge- 
neous Poisson point process on [0, T] with intensity p m G (0, oo), y-a.s. 

Proof. To show this lemma we use Proposition 16.17 of Kallenberg [Kal02j . Ac- 
cording to it, to prove the convergence of H^(X%) to a Poisson point process with 
intensity p m it is sufficient to check that for any interval I C [0, T] 



and 



lim F[H%(X° N ; I) = 0\y] = e~ p ^ 11 (5.26) 



\imsu V E[H%(X N ;I)\y] < p m \I\, (5.27) 

N->oo 

where |/| denotes the Lebesgue measure of /. 

The proof of the first claim is completely similar to the previous ones. We start 
with a one-block estimate for ( 15.261) : 

lim j\rV2j,-i e iV 7 V2/3 2 E[s£(0, U)] = p m , (5.28) 

iV— >oo 

Using the notation of Section [21 we get 

E[^(0,t0] = j A J^y^ Uz ^ ( 5 - 29 ) 

where A m — {z : 3k G {1, . . . , v}Gk(z) > B m }. Dividing the domain of integration 
according to the maximal Gk{z), this is equal 

where D k = {z : G^z) > B m ,Gi(z) < G^{z)\/i ^ k}. Using the substitution 
Zi = hi ± V,iB m on D k (where + sign is used for i < k and — sign for i > k) we get 

NlV ^ mi/ ^yf db l Er=l6f BmGk{b) (531) 
^ (2vr)<V 2 



k=l 



D 
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where D' k = {b : G k {b) > 0, E*=i+i bj + lk-i^Bm > (M ^ k}. The same reasoning 

as before then allows to show that the last expression behaves like p m vN~ x l 2 1'^^^ l 2 ^ 2 
as iV — > oo. 

To compare the real process with the block-independent process, let 

F N (I; X) = l{max{X(z) : iu/r(N) G /} < B m }. (5.32) 

The difference between K[F N (I; X^)\y] and E[Fjv(/; X^)} is again given by the 
Gaussian comparison formula (13.51) . This time the second derivative equals 

5(X(i)-B m )5(X(j)-B m ) 11 l{X(k) < B m } < 5(X(i)-B m )5(X(j)-B m ). (5.33) 

If covariance of X(i) and X(j) equals c, the expectation of the last expression is 
given by the value of the joint density of X(i), X(j) at point (B m , B m ) which is 

(2tt(1 - c 2))-i e -^/(i+c) < c(1 _ c 2)-i exp | _ J 2jV X (5.34) 



(3 2 (l + c) 

The exponential term is the same as in S(c). The polynomial prefactor is however 
different, it diverges faster as c — > 1. We should thus return to (I3.24p with H replaced 
by the right-hand side of f)5.34p . First 

r 1 1 

/ (1-c 2 )- 1 = c _1 argtanh(c) w --log(l-c) (5.35) 
Jo 2 

as c — > 1, which is not bounded for all c as before. The estimates (I3.26P and (13.271) are 
influenced by this change. For (I3.26P we can actually neglect this change, because 
the main contribution to this term came from the neighborhood of d — N/2 (or 
c = 0) and was exponentially small in the neighborhood of d = 1 (or c ~ 1/-/V)- In 
the treatment of (13.271) . the change has more effect, after some computations (13.451) 
becomes 

r 6 ' 

CtN^u- 1 / log(c/ x)e- cNx dx < CtN l/2 v- 1 log N 0. (5.36) 

Jo 

Finally, the change of polynomial prefactor of S implies change in the control of 
(13.281) . The equation (I3.46[) becomes 
</ 

(E2HD < CY,tN~ 3/2 d 2 [l - (l-2dN~ 1 ) 2 p}- 1 exp(Nf(d/N)). (5.37) 

and the linearization of T gives new form of (I3.47P 

CtN*' 2 f xe- c ' Nx dx < CtN- 1 ' 2 0. (5.38) 

Jo 

Therefore, using (15.281) 

pw(4; i) = o\y\ = nF N (i-, x° N )\y] = nF N m x 1 ^ + (i) 

= (1 - E[s$(0, U)]f\ r{N)/v -> e- pm|/| . 

This completes the proof of (15.261) . 

It is easy to check (15.271) . By definition, 

E[H™(X N] I)\y}= Yl ^K(i,X° N )\y}. (5.40) 

i:iv/R<=I 
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Since A?- > for i, j in the same block, E[s#(i, X° N )\y] < E[s%(i, X^)]. Therefore, 

(EH < \I\r(N)/uE[8%(0, U)) = p m \I\. (5.41) 

This completes the proof of Lemma 15.41 □ 

Proof of Theorem I5.M Checking the convergence of finite-dimensional marginals as 
well of condition (i) and the second part of (ii) of Theorem 15. II is analogous as for the 
original clock process S^. We should thus only prove the first part of condition (ii). 
Namely that, for any 77 and e there exist 5 such that 

¥[w §N (S) > V ] <s, (5.42) 

for all N large enough. 
Let 

w f ([r, r + 6]) = sup{min(|/(t 2 ) - f(t)\, \f(t) - f(h)\) : r <t x <t <t 2 < r + 6}. 

(5.43) 

Fix m such that F[S^(T) > rj/2] < e/2, which is possible according to Lemma [5.31 
If H%(X°;[t,t + 6\) < lthen 

^([r, r + 5])< S%(t + 5)- S%(t) < S™(T). (5.44) 

Hence, 

P^flV, r + 6])> r)\iS%{T) < V /2] < n^(X° N] [t, t + 6]) > 2] < Cp m 5 2 . (5.45) 
We can now show (I5.42p . Estimate 

w §N (5) < max{ws N ([r, r + 25}) : < r < T, r = kS, k <E N} (5.46) 

yields 

T8- 1 

nw §N (S) >v\y}<J2 P K([^. (* + 2)5]) > e\y] 



k=0 

T6- 1 

< F[S%(T) > rj/2] + (k + 2)5]) > 2] 



(5.47) 



k=0 

< e/2 + CT<TVm<5 2 < e 
if 5 is chosen small enough. This completes the proof. □ 

Proof of Theorem \l. 6 A Let TZn be the range of the coarse grained process Sjy. Ob- 
viously, for any 1 > e > 0, 

A £ N (t, s) D {TZn fl (t, s) = 0}, (5.48) 

because if the above intersection is empty, then makes less than v steps in time 
interval [te lN , se lN ], and thus the overlap of ajy{te lN ) and a^(se lN ) is 0(u/N). 

If 7?.7v H (t, s) 7^ 0, than there exist u such that S^{u) G (£, s). Moreover, it follows 
from Theorem 15.21 that for any 5 there exist f] such than 

F[S N (u + ri)e(s,t)]>l-5. (5.49) 

This however means that the process o~n make at least rjr(N) steps between times t 
and s and thus the overlap between ajq{te' )N ) and a^(se lN ) is with high probability 
close to 0. 
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Hence P[A^(t, s)\y] is very well approximated by ¥[R N H (t, s) = ®\y). Since 
stable subordinator does not hit points, that is P[3u : Vy/pa(u) — t] — 0, and Sn 
converge in Ji-topology, 

F[n N n (m) = ®\y] P[{v 7//32 ( M ) : u > o}n (s,t) = 0], (5.50) 

which, as follows from the arc-sine law for stable subordinators, is given by the 
formula (Q3I1 . □ 
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